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1.0 SUMMARY OF WORK ACCOMPLISHED 


A computer program has been developed for calculating the radial dis- 
placement due to the body tide as a function of spatial position and time. 
The positions of the moon and sun are pvpI ated by means of the Hill-Brown 
and Newcomb theories, respectively. 

Utilizing Green's functions tabulated by Farrell for the elastic earth 
responses due to point loads, a convolution program has been developed to 
evaluate the elastic earth responses due to ocean loading by a given ocean 
tide constituent. These responses Include radial and horizontal displace- 
ment, tilt and gravity acceleration, and strain tensor components. 

Using the convolution programs, the global radial displacement due to 
the M 2 tide has been obtained. 
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2.0 WORK COMPLETED DURING REPORTING PERIOD 


The tide of the solid earth is composed of two basic responses; (1) a 
body tide due to the yielding of the earth to direct forces of the sun and 
moon and (2) a "load" tide produced by surface forces from the time varying 
ocean tides. It is difficult to distinguish between these two responses 
because their time dependence is similar, being driven by the same ultimate 
force. However, the nature of the driving force of the "body tide" is well 
understood, while the knowledge of the deep ocean tides through global 
numerical modeling is a recent advancement. 

The body tide varies in a relatively smooth nature over the earth's sur- 
face, depending principally on averaged overall elastic properties while 
the load tide is complicated by discontinuities of the surface load at 
coastal boundaries and by local ocean tide circulation (e.g., amphidromes 
and anti -amphidromes). Moreover, the displacement of the load tide is ap- 
preciable only in the crust and upper mantle, while the body tide has 
relatively large amplitude through most of the earth's interior. The load 
tide response then depends more on local crustal properties so that varia- 
tions in near surface earth structure will be more reflected in the load 
tide. 


The work performed during the first six months consisted of implement- 
ing an analytic program for representing the body tide deformation, and 
establishing the required software for evaluating the responses of the earth 
to the ocean loading tide. These tasks will be discussed separately in 
the following sections. 

2.1 Body Tide Program 

The body tide elevation is given by 

M.R ® /Ro\n+l 

Ug {(^, x; t) = I 

e n-t ' a' 
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where 1s the mass of the disturbing body (Moon or Sun), Is the geo- 
centric distance to the body and h^ are Love numbers. The angle 
denotes the geocentric zenith angle of the moon (sun) from the point of eleva- 
tion. The terms In the expansion fall off rapidly so only the first term 
Is of major significance. 

A computer program has been developed to evaluate the solid earth tide 
(body tide) at a user specified position on the earth for a desired time 
Interval; the algorithm presently evaluates for a one day period at one 
hour Increments. Required Input are the Modified Julian Day, and the geo- 
detic coordinates of the point of Interest. 


The radial body tide, Ug , Is evaluated as 



where r represents the unit radius vector at the point of Interest on the 
earth: 


r = \,ao6^ao8\t oos^sinXt ein^'] 

where geodetic latitude, ({> , and longitude, \ , are program Inputs. 

represents the unit vector from the center of the earth In the 
direction of the disturbing body: 

Rj, = [X',u',v'] 

where X', y', v' give the position In earth-fixed coordinates. 

In calculating R for the moon, a true longitude and the latitude 
(above the plane of the ecliptic) are derived from the Hill -Brown theory 
using the Modified Julian Day. Brown's tables express the coordinates of 
the moon as sums of periodic terms whose arguments are algebraic sums of the 
multiples of l, £', F, D, r . See Tables I and II. 
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TABLE I 

Ecliptic Elements 


MJD ■ Modified Julian Day 
D - MJD-241 5020.0 
01 = 0*1. E-4 

z • 296. 104608+1 3. 0649924465*0+0. 0006889*( 01)^ 
£' - 358. 475845+0. 9856002670*0-0. 00001 12*(01)^ 
F = 11. 250889+13. 2293504490*0-0. 0002407*(01)^ 

0 • 350. 737486+12. 1907491914*0-0. 0001076*(01)^ 
r « 281. 220833+0. 470684*01+0. 339E-4*(01)^ 
e * 23. 452294-0. 0035626*01-0. 123E-6*(01)^ 

6g » 99.6904833+360.98564733*0-180.0 


4 


Table II Development of Lunar Position 


Coeff of sine In 6X„ 

n 

(Seconds of arc) 


Multiples of 
£ t' F Dr 


22639.5 

- 4586.426 

2369.902 

769.016 

- 668.111 

- 411.608 

- 211.656 

- 205.962 

- 125.154 

191.953 

- 165.145 

147.693 

- 109.667 


Coeff In sine in Latitude, iji 

n 

(Seconds of arc) 

13461.48 

1010.180 

- 999.695 

- 623.658 

117.262 

199.485 

- 166.577 

61.913 


1 0 
0 0 


1 1 


1 1 
1 1 
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Table III Development of ain and ooa 


Coef. X 10® of cosine in ooa X^ 
and of sine in ain x^ 


Multiples of 
A F Dr 


99972 

1674 

32 

1 

2 

-1675 

-4 

-2 

4 

-4 


0 10 0 1 
0 2 0 0 1 

0 3 0 0 1 

0 4 0 0 1 

0 10 11 

0 0 0 0 1 

0-1 0 0 1 
0 1 0-11 

00 1-10 
0 2-1 12 
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Time (Hours) 



The derived ecliptic lunar position converted to Inertial 


coordinates 


X' ■ aoaXaosd 

m m 


u' " 8in\jsoB^j30Bt - BinBBinc 
m m m 


v' ■ Bin^jsoBt + BinXjooB^-Binc 

m m ^m 


where c Is the obliquity to the ecliptic and 

X ■ i- + D + r ♦ 6X„ 

The conversion to earth-fixed coordinates Is accomplished by a matrix 
transformation 


IL) -(MR.) 

' ^'EF ' Inertial 


where 


oobQ^ Bin&„ 0 
9 9 


M ■ I -ainQ^ cobQ^ 0 

y 9 9 

\ 0 0 1 


and 0 Is the Greenwich hour angle (Table I). 

The unit vector In the direction of the sun Is derived from 
Newcomb's theory In the same manner. The ecliptic elements, t, i', F, D, r , 
are the same as for the lunar development. The solar coordinates 
(ooBX^tBinX^) are expressed as algebraic sums as listed In Table III. 


The conversion to equatorial coordinates Is 


X' ■ 00&X. 


»j' ■ §in\^oc9t 


«inX^«inc 


and a final transformation to earth-fixed coordinates Is applied* using 
Equation (1). 

A comparison of the radial body tide displacement computed by this 
analytic procedure and the ephemerls calculation from GEODYN Is displayed In 
Figure 1. 

2.2 Earth Response to Ocean Loadin 


The calculation of the earth deformation due to surface mass loads 

closely follows the Green's function approach of Longman and Farrell. 

Farrell has Integrated the equations of motion for a self-gravitating 

elastic spherical earth using a Gutenberg-Bullen A earth model and produced 

load Love numbers h', 1', and k' to high order n . The elastic earth 

n n n 

response then reduces to a convolution of the ocean tide with the Green's 
function over the global ocean. The appropriate Green's functions for the 
augmented potential and surface vertical and horizontal displacements at an 
angular distance y (spherical earth) from a point load at the pole per unit 
of loading mass are 


«'(y) “ I 

”e n-0 " " 


U'(y) ■ 5" I 

% n-0 " " 

D • aP^(cc5eX) 
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whtre Is the mass of the earth, R Is the mean earth radius and g Is 
the acceleration of gravity at the surface. Green's functions for the de- 
flectluM of the local vertical and the differential gravity acceleration are 


T'(y) 

G'(y) 


1 Z 3P (o<?*y) 

^ n-0 " " 


t?- I (n+2h;-(n+l)li-) P JooeY) 
% n-0 " " " 


Similar expressions can be obtained for the non-zero elements of the strain 
tensor at the earth's surface. As pointed out by Farrell, these Green's 
functions are slowly convergent series and must be sunned to large values of 
n . 


The sunned values for these Green's functions as a function of angle 
have been obtained from table values given by Farrell, where values are 
available for U', V', T', G' and , where S' Is the strain tensor. 
The other diagonal components of the strain tensor (off-diagonals are zero 
for the symmetric point load coordinate system, with load at the pole, for 
which the Green's functions are derived) are calculated from 



U' * V' 

- ecte — 

m O 



where x(a) and o(a) are the Lame parameters at the top layer of the 
earth model. Here primes on the Green's functions denote that they are In 
the symmetric point load coordinate system. 


The earth responses to ocean tides are then 

• jjp^ C(«'.x';t) pG^(Y)dP' 
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where 


oo$y ■ eineeine' + oo§iao9^'oos{\“X') 

end c(etXit) denotes the ocean tide. Here G^(y) represents the Green's 
function discussed above transformed to an earth-fixed coordinate system and 
R represents the appropriate response. The global ocean tide modeU 
developed for GSFC and available In the CALTOR program are those used In the 
present work. The total tide Is approximated as 

t(a.x;t) ■ y (e.x;t) 


where 


C^(«.x;t) ■ ooa[o^(t-t^)+e^(e,x)+c^] 

and 1 ■ M^t Sg. K 2 » response Is then 

^(♦.x;t) • y fl^(»,x;t) 


where 


i?^(^.x;t) 



soaa^{t~t^) 


0, 


and 


^1g * jj •^^0^<(f'*^')‘»®«C'J'i(^'tX')^c^]G^(Y)«in^'d^'dX' 


1q 


n 


To evaluate the convolutions, It must be realized that the primed Green's 
functions described earlier are with respect to the symmetric point load 
coordinate system, and to resolve components of vector and tensor quantities 
appropriate transformations must be applied. Let (41, x) be the latitude 
and longitude of the point of evaluation for the convolution and ((|i',X') be 
the latitude and longitude of the water column being considered as the load. 
Moreover, let r, 0, x denote the unit vectors In the direction of Increas- 
ing r, 0, X at the point of evaluation (where 0 Is the co-latitude). 

A ^ 

Then the Green's function for the horizontal displacement In the 0 and x 
directions, respectively, are 


Ve(Y) = V'(Y)coao 
Vj^(y) = V'(Y)eina 


where 


0080 . 


sin^'-sin^coBy 

oos^siny 


sino 


ao8(t)'8in{\-A ) 
8iny 


and as before 


008y = + qo 8 ^oo 8 ^'cob{\-\') 


einy = /I -ooB^y 

Similar expressions for the radial displacement and strain tensor components 


are 


U(y) - U'(y) 


S - S' 
rr rr 


2 2 

S ■ 008 a S' + ein a S' 
60 YY 


S.„ ■ Binaaoetx (S' -S ' ) 
X0 ' YY XX 


2 2 

S.. ■ sin a S' + aoe a S.'. 

XX YY XX 


Software is presently complete for evaluating the convolutions for the 
following responses: 


radial displacement, ((|>,x;t)] 


eastward displacement, i?[V (i|',x;t)] 

s- 


southward displacement, /?[V. ((|>,X;t)] 

®i 


rr Strain Tensor component, /?[S^^ ( 4 i,x;t)] 


66 Strain Tensor component, i?[S.. ( 4 ',X;t)] 

00 ^ 


X0 Strain Tensor component, /?[S ((>,X;t)] 

A 0 ^ 


XX Strain Tensor component, i?[S (i}>,x;t)] 

A A ^ 


where i, = M 2 , $ 2 , N^, K 2 , , 0^ , . The output of the software is an 

amplitude, A^(i(‘,x) , and phase, 'i'^((j),x) , defined at each grid point 

over the surface of the globe. The response is then 
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CJOsfo^ (^»X) + E^ ] 

The amplitude and phase have been computed for the global radial displace 
ment due to M 2 and are displayed In Figure 2. 

2.3 Conformance to Schedule 


Work Is currently behind projected schedule. This Is due primarily 
to the slow and unpredictable turn-around on the GSFC 360-91 using STAND-BY 
computer time. Program development work Is done on the 360-95 when practic- 
able, but the 360-91 Is more efficient for the high-core requirements of the 
programs for execution. At this point In the work, we had anticipated 
evaluation of responses for the radial and tangential displacements, tilt, 
gravity acceleration and strain tensor components due to the M 2 tide. 
Currently we have only the radial displacements for the M 2 . We antici- 
pate better turn-around during the holiday season and some of the delay in 
schedule can be made up. 

2.4 Analysis of Progress 

With the exception of delayed turn-around detailed in Section 2.3, 
progress was satisfactory during the reporting period. The computer pro- 
grams were developed and checked-out, and the body tide radial displacement 
subroutines are suitable for stand-alone use. The amplitude and phase chart 
for the radial response to the M 2 tide loading is also of interest as a 
partial result. 

2.5 Efforts to Achieve Reliability 


Not applicable. 

2.6 Personnel Changes 

There have been no personnel changes during the reporting period. 
Persons involved in the work to date are 
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Ronald Estes - Principal Investigator 
Daniel Chin - Senior Programmer/Analyst 
01m Strayer - Programmer 


2.7 New Technology 

There has been no new technology developed, as we understand the mean 
Ing of the new technology clause, during the reporting period. 



3.0 WORK PLANNED FOR NEXT REPORTING PERIOD 

(A) The convolution evaluation program will be run for the M2. $2. K^ 

constituents to produce amplitudes and phases on a global grid for: 

(1) radial displacements (M2 already complete) 

(2) horizontal displacements (east and south components) 

( 3 ) strain tensor components 

( 4 ) tilt and gravity acceleration. 

(B) The amplitudes and phases for the earth responses will be expanded Into 
a series of spherical harmonics In the form 

A ao 8 \b = y [a coemx+b einmx] P _(cosd») 

» m nm nm nm 

n.m 

A sirup = y [c „cosmx+d„ sinmx] P^„(c?c?S(|i) 

*■ nm nm nm 

n.m 

A special purpose program will be fabricated to perform this calculation. 

(C) A software system will be fabricated to evaluate the radial displacement 
due to ocean loading and the body tide as a function of ifi.x and t . 

Effects due to M2. $2. and will be included. The calculations will 

be in terms of spherical harmonic coefficients. 

(D) The software system described in (C) will be designed such that the user 
may substitute the spherical harmonic expansion coefficients for any of 

the earth responses in place of those for the radial displacement and obtain 
the response as a function of 41.x and t . 



Thousands of Dollars 



C. Funds Expenditure Report (Semi-Annual) 


Total Task 

26.8 


Prior Expenditure 


0.0 

Expenditure this Period 


15.2 

Total Expenditure to Date 


15.2 

Amount Remaining 

11.6 


Total Hours Expended 

824 


Amount Next Period 

11.6 


Amount Require to Complete 
Work 

11.6 



